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Abstract 



^ I We calculate the magnetic moments of heavy baryons with a single heavy quark in the bound- 

ed ■ state approach. In this approach the heavy baryons is considered as a heavy meson bound in 

the field of a light baryon. The light baryon field is represented as a soliton excitation of the 
light pseudoscalar and vector meson fields. For these calculations we adopt a model that is both 
chirally invariant and consistent with the heavy quark spin symmetry. We gauge the model action 
with respect to photon field in order to extract the electromagnetic current operator and obtain 
^ , the magnetic moments by computing pertinent matrix elements of this operator between the 

■~«.^ ' bound state wavefunctions. We compare our predictions for the magnetic moments with results 

i-C ' of alternative approaches for the description of heavy baryon properties. 



1 Introduction 



There has been recent interest to study properties of baryons that contain heavy quarks. In particular 
the development of the heavy quark or Isgur-Wise symmetry P has generated many studies on 
^ ' the spectrum of baryons with a single heavy quark in the so-called bound state approach. To set the 

notation for what we consider a baryon with a single heavy quark we denote heavy quarks (charm or 
bottom) by "Q" and light quarks (up or down) by "g" . In the quark language a baryon whose properties 
we wish to explore has the structure qqQ- This defines the quantum numbers of the considered baryon. 
In general, of course, the structure of such a baryon may be more complicated because additional 
quark-antiquark-pairs may be excited. In the above mentioned bound state approach the heavy 
baryon is considered as a heavy spin multiplet of mesons (with structure Qq) bound in the field of the 
nucleon, or more generally a light baryon of qqq-structme. As suggested by the expansion jHUH] 

of QCD the field of the nucleon emerges as a soliton configuration of light meson fields of qq structure. 
Reviews of the soliton approach to the light baryons are given in refs. Cni~C2 while the bound state 
approach for the "light" hyperons is discussed in refs. |15j-|18j. 

The basic idea is to consider an effective meson model that includes both, light degrees of freedom 
and mesons that contain a heavy quark. In such a model the light components of the heavy meson 
fields couple the light meson fields according to the rules of chiral symmetry. Thus this it permits not 
only an expansion of heavy baryon properties in powers of 1/M, 1/Nc |19| but also in the number 
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of derivatives acting on the light components of the heavy system. However, this introduces a large 
number of unknown parameters in the model Lagrangian. We therefore find it more compelling to 
employ a model with light vector mesons (p and uj) to construct the soliton rather than a model with 
many derivatives of the light pseudoscalar fields. Actually the need for light vector mesons is not 
surprising since, in the soliton approach, they are necessary to explain, for example, the neutron- 
proton mass difference [201, the nucleon axial singlet matrix element [2] and the "high-energy" 
behavior of phase shifts in meson-baryon scattering |12j . Furthermore heavy quark components of 
the heavy meson field are required to exhibit the heavy quark symmetry as the masses tend to 
infinity, i.e. spin and fiavor independent interactions. Therefore the model also requires both heavy 
pseudoscalar and heavy vector meson fields. Based on a model |^ that reflects all those features, 
the heavy baryon mass splittings have been discussed |23l . obtaining satisfactory agreement with 
experiment. Also pentaquarks of structure qqqQq have been considered in this model as antimesons 
(Qq) bound in the field of the soliton. 

Here we will go one step further and study magnetic moments of baryons with a single heavy 
quark in the model suggested in ref. j22j as an example for computing static properties. This not only 
requires to construct the soliton of the light meson fields, to compute the heavy meson bound state 
wavefunction, and to quantize the the combined system to generate states to be identified as heavy 
baryons but in addition also to derive of the electromagnetic current 19.12. operator in the model 
and to subsequently evaluate matrix elements thereof between heavy baryon states. It is important 
to note that both the light and the heavy sector of the model will contribute to that current and 
hence to the magnetic moments of heavy baryons. Stated otherwise, these magnetic moments directly 
vary with the parameters of the light sector and not only indirectly via the soliton profiles as e.g. the 
spectrum of the heavy baryons does. 

This paper is organized as follows. In section 2 we will review the model Lagrangian for both 
the light and heavy flavor sectors of the model. This will also provide the opportunity to derive 
covariant expressions for the electromagnetic current. Section 2 also contains a brief discussion the 
soliton profiles that emerge in the light flavor sector as well the heavy meson bound state profiles. 
In section 3 we will generate states with good baryon quantum numbers by canonically quantizing 
collective coordinates that parameterize the (zero-mode) rotations of the combined soliton - heavy 
meson system. This makes mandatory the discussion of field components that are absent classically 
but get induced by the collective rotation. Section 4 represents the major progress for the studies of 
heavy baryons as a heavy meson bound to a soliton as we construct the magnetic moment operator 
in this model and compute the corresponding matrix elements numerically. Concluding remarks are 
to be found in section 5. Also we complete the paper by including an appendix that summarizes 
functionals of the meson profiles that emerge along the computation. 

Some preliminary results of this study have already been presented in a conference contribution |24| . 

2 The model Lagrangian 

The model Lagrangian can be divided in two distinct parts. The first part, >Ciight describes the 
interaction of the mesons that are built out of the light quarks up and down. Our model not only 
contains the pseudoscalar pion fields but also the vector mesons co and p. The interactions among 
these light mesons is dictated by chiral symmetry. The second part, £hcavy hi addition contains fields 
that correspond to mesons that are built out of a single heavy quark (charm or bottom) and light 
quarks^. This part of the Lagrangian is consistent with the heavy quark spin flavor symmetry and 
therefore contains both pseudoscalar and vector degrees of freedom. In the charm quark sector these 
are L»(1865) and i:>*(2007) while 5(5279) and B*(5325) contain a bottom quark. 

^In what follows we will refer to mesons with a single heavy quark as heavy mesons. 
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The Lagrangian Cnght contains (topological) soliton solutions that are identified with baryons 
according to the Skyrme-model picture. Substituting this soliton configuration for the light meson 
fields in £hcavy generates a potential for the heavy meson fields. The Klein-Gordon equation associated 
with this potential has solutions with energy less than the rest mass of the considered heavy mesons, 
i.e. a bound heavy meson. We combine bound heavy mesons with the soliton of the light mesons (i.e. 
the soliton represents a baryon of light quarks) to obtain representatives for baryons with a single 
heavy quark. In the soliton description this model for the heavy baryon is that of a heavy meson 
bound to a light baryon. In this section we will review how this picture emerges and also discuss the 
basics for computing the magnetic moments of such heavy baryons. 



2.1 The light meson Lagrangian 

The light sector of the model is given by the chirally invariant Lagrangian discussed in refs. |25ll28j . In 
addition to the standard Skyrme model, which considers pseudoscalar pion fields only, this Lagrangian 
also contains vector meson degrees of freedom. In order to efficiently incorporate chiral symmetry it is 
most useful to employ the non-linear representation of the isovector pion fields (tt) via the chiral field 
U = exp {in ■ t/Ztt) where = 93MeV and f are the pion decay constant and the Pauli matrices, 
respectively. To begin with, a chirally invariant Lagrangian for vector mesons should contain both, 
vector and axialvector fields. The chirally covariant elimination of the axialvector introduces the root, 
£, of the chiral field, i.e. ^"^ = 11 |2H| as well as vector and axialvector currents of the pion fields: 

V, = \ {id^^^ + ^^5^0 , V^ = \ {^d.i^ - i^d.i) . (1) 

We comprise the isoscalar lo and the isovector p within a single matrix field ~ (p^ •T + a;^]|)/2 and 
further define ^ 

Ft_,i,{p) ^ df_,py - dyp^ + ig[p^,py] and i?^ = p^ - -w^^ (2) 

to compactly write the normal parity part of the Lagrangian for the light mesons 

Cs = fllT [p^p^] + -mlfliT [C/ + C/t - 2] - itr [FMF^^^p)] + m^lT [R,,R^^] . (3) 

Here — 138MeV and my — 770MeV are pion and vector meson masses, respectively. The last 
term in eq. iPJ contains the tri~linear vertex p^ ■ (7? x d^i:) and allows to determine the coupling 
constant g « 5.6 from the decay p — *■ tttt. In addition the light meson Lagrangian contains a part 
that involves the Levi-Cevita tenor e^^^^. Since this part includes the nonlocal Wess-Zumino-Witten 
term it is most convenient to write it in action form using differential forms like p = p^j,dx^ 



r,.,{p,R,F) = / trb'] 



tr 



^71 + ^72) V - ^l2F{p){pR - Rp) - 2tg\j2 + 2j3)R^p 



Note that we will only consider the physical case of Nc = 3 color degrees of freedom. 

In ref. |25j two of the three unknown constants 71,2,3 were determined from purely strong interaction 
processes. Defining h — —^^71, gvvct> = 912 and = ^ the central values h = 0.4 and gvvct> = 1-9 
were found. Within experimental uncertainties (stemming from the uncertainty in the uj — <j) mixing 
angle) these may vary in the range h ~ —0.15, .., 0.7 and gvviti ~ 1-3, •., 2.2 subject to the constraint 
\gvvit> — h\ K 1.5. The third parameter, k could not be fixed in the meson sector, however, from the 
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study [2ni of nucleon properties in tlie U{2) reduction of the model it was argued that k w 1. To be 
specific, we will always employ the parameters: 



g = 5.6, h = 0.3, gvv0 = 1- 



K = 1.2. 



that have been found suitable for describing the low-lying baryons for three light flavors |3(J| . 
The total action for the light mesons is the sum 



(4) 



(5) 



The electromagnetic current associated with Fught is most easily obtained by gauging it with the 
photon field Af^ Q where Q — diag(2/3, — 1/3) is the quark charge matrix and extracting the linear 
term: 

flight ^ ST[^,p,A 



This yields 



(6) 



J 



light 



f> { Q [e^p.e - ep.el } - ^tr { Q [^R,e + eti?^e] } 
I ^ 3 " 



2i 




3" 


.4772 " 






1 


71 + 


9 





71 



:72 



wtr{[e^Qe + eQe^]pVp"} 



72 e^.p^tr { [C^ - CQC^] [^»" + P V^?" - p'R^pn } 
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F^PR" + R^PP" - - [F^Pp" + p'^FP" 
9 

F^PR" + R^FP" + - (F^^pP + p^FP" 
9 

R'^RPR" - - [R^RPp^ +p''RPR'' - R'^pPR"] 
9 



-ig' (72 + 273) e^,p,tr <j Q ( C 
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R'^RPR" + -{R^RPp" + p'^RPR" - R^pPR" 
9 

+2die^,,,a''{tr {Q " p'R") + - P^i?") } • 



(7) 



The last term, proportional to di has no analogue in the pure hadronic part of the action (@J . This 
term is of electromagnetic nature and the coupling constant can be related to the light vector meson 
decay uj tt^: \2di - ^\ = 0.038 ± 0.002 



2.2 Solitons 

The action for the light degrees of freedom contains static solition solutions. The hedgehog type 
configuration reads jSH): 



U{f) 
Pl{r) 



exp {if ■ f F (r)) ,i.e 



^(7^) — exp [if-f 



gr 



and i^{f) 
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,F{r 



(8) 



while all other field components vanish classically. These configurations are invariant under combined 
spatial and fiavor rotations generated by the vector sum L + S + 1, the so-called grand spin. Assuming 
the boundary conditions (a prime indicates the derivative with respect to r) 



f(0)=7r, G(0) 



w'(0) = 0, F(oo) = G'(oo) = cj(oo) = 



(9) 
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yields baryon number one solitons by substituting the ansdtze (jHJ into the action @ and extremizing 
the resulting functional of the classical mass, Mc\[F, G,uj\. For completeness we list the explicit form 
of that functional in the appendix. Note, that this soliton does not yet describe baryon states with 
good spin and flavor quantum numbers. Before describing the cranking procedure that generates such 
states, we will focus on that part of the model Lagrangian that contains the heavy meson fields. 

2.3 The relativistic Lagrangian for the heavy mesons 

The relativistic Lagrangian iZheavy which describes the coupling between light and heavy mesons is 
given by |22j : 

+ 2iMd {Pp^^Q^^ - QmP^'^O - ^e"'"'' [Q^^c.P^.Ql + QpP^iQ.a)^ 
- ^^^^ |2Q^F^'^(p)Qt _ [D0PF,Ap)Qi + Q^F,Ap){D^P)^ | . (10) 

Here the mass M of the heavy pseudoscalar Meson P differs from the mass M* of the heavy vector 
meson Q^. We take 

M = 1865 MeV (D - meson) , 5279 Me V (B - meson), 
M* = 2007 AfeV (D - meson) , 5325 MeF (B - meson) . (11) 

It should be noticed that the heavy meson fields are conventionally defined as row vectors in isospin 
space. The chirally covariant derivatives of the heavy fields P and Q and the field-strength tensor 
Q^i, are defined as: 

D^,P^ = [d^ - iagp^ - i {1 - a) Vf_,] P^ = [d^ - iv^ - igaR^] P^ , 

D,Ql = [d^-iv,^igaR^]Ql, (Q^,)^ = D^Ql - D^Ql (12) 

In eq. H10|l we have two sets of two terms each that involve the identical coupling constants d and c. 
Each of these terms is chirally invariant and could thus in principle carry independent coupling con- 
stants. However, the condition that the Lagrangian £hcavy obeys the heavy quark spin symmetry in 
the limit M, M* — > oo requires the combination of these terms as given in eq. (Uni 22. Note however, 
that we do not assume infinitely large masses for the heavy meson. The coupling constants d, c and 
a are still not precisely determined. While d and c can be determined from the decay widths of the 
heavy mesons 

d = 0.53, c=1.60, (13) 

there is no direct experimental information for the value of a. The value of a would be unity if a 
possible definition of light vector meson dominance for the electromagnetic form factors of the heavy 
meson was adopted [221 ■ 

Due to the electromagnetic gauging of the Lagrangian extra terms appear in the light vector fields 
Pfi, Vfi and Pf^. The respective electromagnetic and chirally invariant forms read (we choose the electric 
charge to be unity): 

Pm=Pm--2A<- v^.^v^ + ^A^,{^Q^^ +^^Q^-2Q) and = p,, - (^2^^ - C^Q?) . 

_ ^ (14) 

Considering these extra terms, the electromagnetic and chiral covariant derivative of the heavy meson 
fields has to be modified: 

D^P^ = {d^-iagp^-i{l-a)Vf,-i{Q-C)A^)P^ 
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= - iagpfj, - i {1 - a)v^ + ieA^ 

which makes plausible the above assertion that a = 1 corresponds to vector meson dominance because 
for that value there is only the direct coupling to the heavy mesons. The charge of the heavy quark in 
question is denoted by C, ie. C = | for the charm sector and C ~ for the bottom sector. Finally 
we substitute eq. (fTB|l . for D^P\ its analogue for [Df^Q^Y as well as p, etc. from eq. (fnji for p into 
the Lagrangian (|10|l to obtain the electromagnetic current of the heavy mesons 

= i (pC (D^P)^ - D^PCpA + I (q^CQI, ~ Qt^uCQ'^A 

^ =0 

-leMd [P [m^ - QO Ql ~ Q, {iQS} - ^ QC) P) 

(pCFP^{p)Q^^ - Q^fP^{p)CpA . (16) 
my \ I 

Here the abbreviation C = C — (■^Q'^^ +C^Qi^) has been employed. A similar electromagnetic 
current was also derived in ref. p9] to discuss radiative decays of heavy vector mesons. However, the 
d-type term in eq. H16(l was omitted as it does not contribute to such processes. 



(15) 



/T- heavy 



2.4 Bound states 

Substituting the soliton configuration for the light meson fields in /Zhoavy generates a potential 
for the heavy meson fields. The existence of heavy meson bound states in this potential has been 
established some time ago j23) . Such bound states emerge for both, positive and negative frequency 
modes. The former correspond to mesons bound to the light quark soliton and thus carry quantum 
numbers of ordinary three-quark baryons. The latter, however, have antimesons bound to the soliton 
and their quantum numbers cannot be constructed in a three-quark picture. Rather additional quark- 
antiquark excitations are required and the resulting heavy baryons possess pentaquark structures. 

In the current study we want to compute the magnetic moments of ground state baryons with a 
single heavy quark, i.e. Ac,Sc when the heavy meson is a D-meson and Ab,I]f, in the case of the 
B-meson. The ground states, of course, correspond to the most tightly bound state. Due to the 
hedgehog structure of the background soliton, which has nonzero orbital angular momentum, these 
bound states emerge in the P-wave channels^, rather than in the S'-wave as nai'vly assumed. For that 
reason we will only consider the P-wave channel and refer to ref. |23| for details on S'-wave channel 
bound states. The appropriate ansatz for the heavy meson multiplet in the P-wave channel reads: 



Pt(r,t) = 



1 



47r 
1 



$(r)f.Tp(e)e"*, Ql{f,t) 



1 



p(e)e-*. 



^o{r)pie) e*' 



(17) 



Here p(e) 



P2(e) 



refers to a properly normalized isospinor describing the isospin of the the heavy 



meson multiplet. Upon canonical quantization the Fourier amplitudes p(e) and p^(e) are elevated to 
annihilation and creation operators for a heavy meson quantum with the energy eigenvalue e (see 

■^The orbital angular momentum quantum numbers denote those of the pseudoscalar component Pt of the heavy 
meson multiplet {Pt,Q]l). 
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below). We substitute the field configurations ^ and If 7|) into the Lagrangian given by eqs. and 
(|f 0|) and integrate over coordinates to obtain the Lagrange function of the form 

L = -M,i[F, G, Lu] + h[F, G, u; $, ^-o, ^i,^2W ie)p{e) . (f 8) 

Here Md denotes the classical solition mass. Its variation yields the profile functions F, G and uj. The 
heavy meson fields are contained in the functional where the subscript indicates the parameterical 
dependence on the frequency of the fluctuation heavy meson fields. The explicit expressions for the 
functional Md and are given in the appendix. We vary with respect to $, ^o, ^"1 and ^'2 to get 
the equations of motions for the heavy meson fields. These constitute Klein-Gordon type equations 
with potentials generated by the profile functions F, G and uj. We then tune the frequency e such 
that these equations yield a normalizable solution with |e| < M. The solution with the smallest 
such |e| is the bound state we are looking for. The equations of motions for the heavy meson fields 
are homogeneous linear differential equations. Hence they provide a solution only up to an overall 
prefactor. Nonetheless, the equations of motion for the heavy meson fields allow us to extract a metric 
for a scalar product between different bound states. In particular its diagonal elements can be used to 
properly normalize the bound state wave functions. The physical interpretation of the normalization 
condition is that each occupation of the bound state should add the amount |e| to the total energy 
and that such a bound states carries unit heavy flavor. Thus we obtain the normalization condition 



1 (19) 



in addition to the canonical commutation relation [pi(e), p](e')] — 5ijS<i,t' for the Fourier amplitude of 
the bound state. For further details we again refer to ref. pS] . 



3 States with Spin and Flavor Quantum Numbers 

It is easy to verify that the field configurations for both the light mesons and the heavy mesons are 
neither eigenfunctions of the spin nor the isospin generators. It is the aim of this section to construct 
such eigenstates for the bound state system. 



3.1 Collective coordinates 



We employ the cranking procedure in order to generate states that correspond to physical baryons. 
In a first step we introduce collective coordinates that parameterize the (iso-) spin orientation of the 
meson configuration, 

e(r ,t) = A{t)^H{r)A^{t), p,{f,t) = A{t)p, H{r)AHt)- (20) 
The time dependence of the collective rotations is measured by the angular velocities fl : 

AHt)f^A{t)='-f-n. (21) 

In addition to the collective rotation of the soliton configuration, classically vanishing field components 
are induced. In the case of the light vector mesons these are given by ,26 : 



(22) 



After introducing these additional fields the Lagrangian for the light mesons now contains a term 
which is quadratic in the angular velocities. The constant of proportionality defines the moment 
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of inertia a^[-F, G, w; ^i, ^2, </']■ Varying this moment of inertia with respect to the radial functions 
C2('') and (p{r) yields their equations of motion that are linear, inhomogeneous second order 
differential equations with the inhomogeneity defined by the classical profile functions G and uj. We 
solve these differential equations subject to boundary conditions that avoid singularities in both the 
respective equations of motion and the moment of inertia |2fi| . The moment of inertia as evaluated 
for this solution is a pure number. For the parameters listed above eq Q we find = 5.01 GeV~^. 
The heavy meson fields also need to be rotated in isospin space, this is done analogously to eq. H2U|) : 

P^^A{t)P\ Ql^A{t)Ql. (23) 

We have now summarized all the collectively rotating fields and are in the position to give the Lagrange 
function for the collective coordinates Lagrangian >Ciight + 'Chcavy, 

2 

icon = -M,z + + + fp^f^ • rp (24) 

that is obtained by integrating flight +'Chcavy with the above ansdtze substituted. The new quantity is x 
which describes the coupling between the collective rotations and the bound state wave functions. For 
reasons that will become obvious below, we will refer to x as the hyperfine parameter. It is a functional 
of all radial functions, i.e. it contains both light and heavy meson fields and can straightforwardly be 
computed once a value for a is chosen, cf. tabled 

The quantization of the Lagrangian (|24|l proceeds along the bound state approach to the Skyrmion 
SH]~SH]: Noether charges for spin and flavor have to be constructed. Considering the variation 
of the fields under infinitesimal symmetry transformations, we find for the isospin transformation 
transformation: 



. ' 2. 



--D,,iA)— + .... (25) 

Here $ refers to any of the given iso-rotating meson fields and the ellipses represent subleading terms 
in 1/Nc, e.g. time derivatives of the angular velocities which might arise from eq. H22(l . Furthermore 
Dij{A) = 1/2 tr[TiATjA^] denotes the adjoint representation of the collective rotations A. It is 
straightforward to conclude from eq. (|25|l that the total isospin is related to the derivative of the 
Lagrange function with respect to the angular velocities: 

The total spin operator J contains the grand spin operator G, 

G, = J, + DT^\A)I, = J,; - , (27) 
with J*°' = dL/dCt. As a consequence of eq. (|26|) we obtain the relation: 

(J'"'')2 =/2 = /(/ + !). (28) 

By construction the light meson fields do not contribute to the grand spin G. Also those parts of the 
heavy meson wave functions (|23(l that are placed between the collective coordinates A and the spinor 
p carry zero grand spin. With the normalization condition (|19|l one therefore finds for the grandspin 
operator: 

G=-p^^p. (29) 

This relation connects the operator multiplying the hyperfine parameter in the collective Lagrangian 
in eq. 124|) to the spin and isospin operators. The collective coordinate piece of the Hamiltonian is 
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then obtained from the Legendre transformation: 

i/eoii = n ■ - Lcou = ^ + xg) ' . (30) 

Since the moment of inertia, is of order Nc, this operator is of order l/Nc- Finally the mass 
formula for baryon B with a single heavy quark reads: 

Mb = M,i + |e| + ^ [xJ{J + 1) + (1 - x)I{I + 1)] • (31) 

Contributions of O(x^) have been omitted for consistency since these terms are quartic in the heavy 
meson wavefunctions and have been excluded form the model form the very beginning. Obviously the 
parameter x characterizes the hyperfine splitting. In ref. j23| it has been verified that it vanishes in 
heavy quark limit M, M* — > cxd . This is a consequence of heavy quark spin symmetry that predicts 
hadrons to be degenerate that only differ by their spin quantum number. 

To constrain the unknown parameter a we compute the mass differences of the bound P — wave 
heavy baryons Sc and to Ap. In addition, we calculate the mass differences of Ac to the nucleon 
and Kb- A more detailed comparison containing also S-Wave bound states and radially excited states 
can be found in ref. |23| . For quite a range of a fair agreement with the empirical data is obtained. 
The Sc — Ac mass differences suggest a negative value of a while the experimental mass difference 
between the nucleon and the Ac is reproduced for a ~ 0.4. The Kb — Ac mass difference is off by 
only about 5%. 

Table 1: Quantities that enter the mass formula H31|l as functions of the unknown parameter a. The 
other model parameters are as in eqs. (^1 lllll and p3(l . The binding energy is defined as t^; = M — |e|, 
M being the mass of the heavy pseudoscalar meson. We also compare the resulting mass differences 
to experimental data in cases the latter are available [311 . The superscripts denote the heavy quark 
flavor under consideration. Energies are measured in MeV and masses are displayed relative to Ac. 



a 


-0.3 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


exp. 




615 


548 


526 


504 


483 


461 


439 


418 
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376 
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773 


748 


724 


699 
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651 


626 
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0.165 


0.140 


0.132 


0.123 


0.114 


0.105 


0.096 


0.087 


0.078 


0.069 






0.060 


0.053 


0.050 


0.046 


0.045 


0.042 


0.039 


0.037 


0.034 


0.031 






167 


171 


173 


175 


177 


178 


180 


182 


184 


186 


168 


AM* 


216 


214 


213 


212 


211 


210 


209 


208 


207 


207 


235 




-1187 


-1252 


-1274 


-1295 


-1316 


-1337 


-1358 


-1379 


-1399 


-1419 


-1344 


AMa^ 


3164 


3172 


3173 


3176 


3178 


3181 


3182 


3185 


3188 


3191 


3339±9 



In the framework of collective coordinate quantization baryon wavefunctions emerge as Wigner-£) 
functions of the collective coordinates. These wavefunctions enter the computation of the magnetic 
moments to be discussed in the next section. For baryons without heavy flavor components these 
wavefunctions read 

*^(/ = J, /3, J, ,h) = Md'^1=:^1^ {A) , (32) 

while for the baryons with a single heavy quark component we need to couple a diquark state (J''"' 
is integer) of collective coordinates with unit occupation of the bound state that carries half-integer 
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spin (J - J'°^ = G) 



*B (/, /3 , J, J3) = AA' ^/ jt' , i G3 ^t-J- (^) 1 1 ' ^3> 



(33) 



to total spin according to eq. H27|l . In the above equations, C is a Clebsch-Gordan coefficient while 
M and jV' are suitable normalization constants. 



4 Electromagnetic current and results for magnetic moments 

In this section we compute the magnetic moments from matrix elements of the electromagnetic current 
operator given in eqs. Q and The defining equation for the magnetic moment operator reads 



1 



flight j^hcavy 



(34) 



We now substitute the field configurations 1©, ((13, and (j^Sl into eqs. and We then 

obtain the operator as sum of terms that are products of integrals over the profile functions and 
operators that act in the space of the collective coordinates {A) and/or as creation and annihilation 
operators (p) for the heavy meson bound state, 



(35) 



We omit the energy argument of the spinor p, it is understood to be the bound state energy e The first 
two terms in eq. (|35|l do not contain the bound state wavefunction and are those already considered 
for the two flavor version of the model in ref. ■ Note that for convenience we factorized the moment 
of inertia, for the isoscalar piece to make contact with the notation in refs. From eq. {7)) 

we find for this isoscalar piece. 



9a2 







1 drr^< 


[ 9^ ^ + 



F' sin^ F 



25 
1 



72 (sinF(G' - Ci) + F'[G{1 - Ci) + 6]) 



+ 7r (72 + 273)^^' (6 - 1 + cos F)(l + G - cos F) 
2.9 



F'i^i + 6) - - sin i^(2 + G - Ci - 2 cos F) 
r 



(36) 



while the explicit form of the isovector piece reads 
SttMn 

/ arr' 

1 



J drr^ I /2 sin^ i^- ^ cos F(l + G-cosF) 



3 1 

(71 + -j2)ujF'sm^ F~ —^2 + 2-f3)F'LucosF(l + G- cosF) 
52 2g 

—72 sin F[w'(l + G - 2 cos F) - G'w] 
2ff 



+2di 



F'uj H — cjsinf cosi^ 
r 



(37) 
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For the parameter set (@J) we find numerically fis.o — 0.43 n.m and ^v.o — 6.70 n.m.. On the other 
hand fxs.i and fiv.i are bilinear in the heavy meson bound state wavefunctions. Again, there are 
isoscalar 



8Mn 



^-0 (i?a^'0 + - ^) '^*2) + *2 ( 1 + ^ ) - *1 {%r + *2 + 



F'r 



V2c 



-'i'o *i sini^ - *2-;^ (*2t^'7- - 2*oG' 



mvg 



(38) 



and isovector contributions 
2Mn 



fJ'V,l 



2 



r^'2 + *2 1 



2 



«'l(«'2r + *2 + -Ra*l) 



%/2c 



+cf^'o*i sinF $ (*2wV - 2*oG' 



j drr^\r (1'o*2 - ^0*2) + l'o'J'2 + 2i?„*o*i 

— j ^i^^a - 2MrsinF*2*| 



-2r\€ 



(39) 



with Ra = cosF — 1 + a(l + G — cosF). The factor A/at — 940MeV enters because we measure 
the magnetic moments in nucleon magnetons [n.m.]. Numerical results for both the charm and the 
bottom sectors are displayed in table El 



Table 2: Numerical results for the coefficients in the magnetic moment operator H35|l that contain 
the bound state wavefunctions. Data are given units of nucleon magnetons (n.m.). See also tabled 



a 


-0.3 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 0.7 


0.8 




Charm Sector 


MS,i 


-0.168 
-0.037 


-0.186 
-0.037 


-0.191 
-0.038 


-0.197 
-0.038 


-0.203 
-0.038 


-0.208 
-0.039 


-0.214 
-0.039 


-0.219 
-0.040 


-0.225 
-0.040 


-0.230 
-0.041 




Bottom Sector 




0.070 
-0.006 


0.063 
-0.006 


0.061 
-0.005 


0.059 
-0.004 


0.057 
-0.003 


0.054 
-0.003 


0.052 
-0.002 


0.050 
-0.002 


0.048 
-0.001 


0.046 
-0.001 



We now sandwich the operator (|35|l between the states listed in eqs. H32(l and (|33|l to obtain the 
magnetic moment of a given baryon as linear combinations of the functionals /J.s,Oi This 
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M = M* [MeV] 



Figure 1: Coefficients ^s,i a-nd in the magnetic moment operator 1)35(1 that are bihnear in the 
heavy meson wavefunction as functions of the heavy meson masses for a = in the charm sector, i.e. 
forC = 2/3. 



yields [Elin]: 

m(S§=0-^Aq) = -i^^o-i^^i. (40) 

The quantity jls.i = XMs.o — Ms,i arises from the quantization rule a'^Q = J^°' + xG. For the 
proton and the neutron only ^5^0 and ^v,o enter. We immediately get = 2.45 n.m. and ^J'{n) = 
— 2.02n.m. which agrees reasonably with the experimental data = 2.79 n.m. (exp.) and /^(n) = 
— 1.91 n.m. (exp.). As is well known [^H], the isoscalar combination ^5^0 — + mC*^) is predicted 
somewhat too small in this model, 0.43 n.m. vs. 0.88 n.m. (exp.) while the isovector combination 
(2/3)/iy,o — ~ is off by only a few percent, 4.47 n.m. vs. 4.70 n.m.(exp.). 

In figure n we show the numerical results for the coefficients fis,i and ^v,i as functions of (equal) 
heavy meson masses. Obviously these coefficients tend to zero as the mass increases'^. This behavior is 
expected from the heavy mass limit since properties, such as the magnetic moments, that are related to 
spin of the heavy quark are suppressed by inverse powers of the heavy quark mass. Note that also the 
hyperfine parameter x tends to zero in the heavy limit |23| and thus the reduction of 113,1 implies that 
of /is,!. To nevertheless obtain a non- vanishing result in the heavy limit the authors of ref. [H2I added 
an extra term with an undetermined parameter to ^Choavy that is linear in the photon field and thus 

•^Note that the displayed ^iv,l is multipHed by a factor 10 for clarity. For Af = M* > 6GeV the value fJ,v,i = is 
consistent with the numerical accuracy. 



1 1 

1 1 

2 1 1. 1 



2 1 1. 1 
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does only contribute to electromagnetic properties of the heavy baryons It was argued in ref. |H2| 
that this term would be necessary to describe the radiative decays Q — > P7. This term introduces an 
additional, so far undetermined parameter which may be interpreted as the intrinsic magnetic moment 
of the heavy meson field. Adopting a canonical value, l/2Af, leads to /iA^ ~ 0.37n.m. [3^, a value 
considerably larger than our result. It is important to note that the model we consider is capable 
of describing such radiative decays [221 because it is formulated for finite heavy meson masses rather 
than in the strict heavy limit. Thus there is no need for any additional photon coupling in our model. 
The decrease of the /^s.i a-nd /iy.i with the increase of the heavy meson mass can also be observed 
from table |2] where we list these quantities for the physical cases as functions of the undetermined 
parameter a. In this case, of course, the pseudoscalar and vector meson masses are different, c/. 
eq. (fTT|l . 

We see that the coefficients /is.i and /iy.i that are bilinear in the heavy meson wavefunctions 
exhibit quite a pronounced dependence on the undetermined parameter a. However, this dependence 
does not propagate to the magnetic moments of the heavy baryons as seen in tables O and 0] One 
reason is that the total magnetic moments are dominated by /ig.o and /iy,Oi the coefficients that contain 
only the light mesons profile functions. The other reason is that the combination fis,i — XMs.o ~ Ms.i 
essentially stays constant due to the decrease of x with a, c/. tabled 

Table 3: Magnetic moments of heavy baryons in the charm sector given in units of nuclear magnetons 
[n.m.] 



a 


Ac 










-0.3 


0.12 


2.46 


0.25 


-1.96 


-2.26 





0.12 


2.45 


0.25 


-1.96 


-2.26 


0.3 


0.13 


2.45 


0.24 


-1.96 


-2.26 


0.6 


0.13 


2.45 


0.24 


-1.96 


-2.26 



Table 4: Magnetic moments of heavy baryons in the bottom sector given in units of nuclear magnetons 
[n.m.] 



a 


As 






^B 


S^^As 


-0.3 


-0.02 


2.52 


0.29 


-1.93 


-2.24 





-0.02 


2.52 


0.29 


-1.93 


-2.24 


0.3 


-0.02 


2.52 


0.29 


-1.94 


-2.24 


0.6 


-0.02 


2.52 


0.29 


-1.94 


-2.23 



In table El we show the results for the magnetic moments obtained in bound state approach to the 
Skyrme model. For the charm sector these results are the same as those denoted'' "SET 11" in table 7 
of ref. 221. The calculation in the bottom sector is that of the kaon sector with the replacement 
Mk 5279MeV. Besides the lower scale for the magnetic moments of the light baryons, the essential 
difference is that the contribution of the heavy meson bound states to the magnetic moments does 
not decrease as the heavy meson mass increases. Of course, that is expected, as the model of ref. |17| 
does not reflect the heavy quark spin symmetry. 

We finally would like to compare our predictions with those of other descriptions for heavy baryons. 
Analyses of QCD spectral sum rules yield (231: = (0.15 ± 0.05)n.m., — (2.1 ± 0.3)n.m., 

■^Note that the authors of ref. 1171 normalize their computed magnetic moments to that of the nucleon. 
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Table 5: Magnetic moments (in nuclear magnetons) of heavy baryons using the bound state approach 
to the Skyrme model 17 . In this model we find /ip — 1.97n.m. and /i„ = — 1.24n.m.. 





Aq 










Q = C 


0.20 


1.96 


0.42 


-1.12 


-1.67 


Q = B 


-0.21 


2.22 


0.56 


-1.11 


-1.54 



fi^+ — (0.23 ± 0.03)n.m., and /is" = (—1.6 ± 0.2)n.m.. These values agree nicely with our results. 
On the other hand, this result for the magnetic moment of Ac is smaller than the one obtained from 
assuming a canonical intrinsic magnetic moment of the heavy meson fields |34| . In our approach 
those four magnetic moments {pA^, ^^'^ Ms") contain the isoscalar piece 11^13=0 which 

is dominated by the light meson isoscalar contribution, fisfi- We therefore compare our results for 
the transition magnetic moments, ^^13=0 , which do not contain 11^13=0, to results from light cone 

QCD sum rules /ze^^Ac = — (1-5 ± 0.4)n.m. and /is^^Ag = — (1.6 ± 0.4)n.m.. We find that our 
predictions are only slightly larger (in magnitude), in addition, those sum rule results have sizable 
error bars. 



5 Summary 

In this study we have employed the bound state approach to compute magnetic moments of heavy 
baryons with a single heavy quark. In this approach the heavy baryon is constructed from a heavy 
meson field that is bound in the background field of a light baryon. The latter emerges as a soliton of 
light meson fields. In the model, that we consider here, this soliton contains both light pseudoscalar 
and light vector meson fields. This extension of the original Skyrme model is known to reasonably 
describe the phenomenology of light baryons. This is particularly the case for the magnetic moments 
of proton and neutron. For the heavy sector we have adopted a relativistic Lagrangian which does 
not directly reflect the heavy quark symmetry. Rather this model embodies the physical values of 
the heavy meson masses. However, in the limit of infinitely large heavy meson masses it properly 
reflects the heavy spin flavor symmetry. We have then introduced and canonically quantized the 
collective coordinates that parameterize the spin-flavor orientation of the soliton and the bound state 
wavefunction to generate baryon states. 

We have extracted the operator of the electromagnetic current from the electromagnetically gauged 
action. To compute baryon magnetic moments we have sandwiched the appropriate combination of this 
operator between baryon states. For the magnetic moments of the heavy baryons experimental data 
do not yet exist, thus we have compared our results with to predictions of other models. Specifically, 
results are available for spectral sum rule and QCD light cone sum rule analyses. Although our 
prediction for the transition magnetic moments, /^„i3=o . is slightly larger than in the QCD sum 

rule approach, the overall picture is that the predictions in these two descriptions for heavy baryon 
properties are very similar. On the other hand, we have seen that the bound state approach to the 
SUf{N) Skyrme model, which does not respect the heavy spin flavor symmetry, gives significantly 
different predictions to the magnetic moments of heavy baryons with a single quark. Not surprisingly, 
these differences are quite distinct for the magnetic moment of Aq which is most sensitive to the heavy 
quark contribution to the magnetic moments. 
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Appendix 

In this appendix we list the explicit expressions for the functionals of various meson profiles that 
parameterize the collective coordinate Lagrangian H24(l . These functionals have already been reported 
in the literature |25l I14L 1^ . However, because notations do vary within those papers, we list them 
here for completeness. 

First we focus on quantities that only involve the light pseudoscalar and light vector meson profiles, 
cf. eq. The classical mass is given by 



Mr: 



An J I Y iF''^r^ + 2 sin F) +mlf^{l- cos F) 
2g^ 2r^ 



+ (1 + G - cos F)2 + IIf'c^ sin2 F-^G'cu sin F 
r 9 9 

+—F'ujG(G + 2) + i (72 + 73) F'uj \l - 2(G + 1) cosF + cos^ F] !> . (41) 
9 9 ' 

Its variation gives the soliton profiles F, G and uj. The moment of inertia reads 



a' = ^ fdr{fysm^F-±i^'' + 2^+mW) + ^AiCi+C2r + 2{^,-l + cosFr 
1 



3 J I 2g 



+ ^[(3ef + 2C,e2 + + 2(G2 + 2G + 2)e| + ^G^iC! + 66 - 26 - 6 + 1)] 

+ -7i(pF' sin^ F + -j3ipF'[{G - 6)(1 - cosF) + (1 - cosF)^ - G6] 
9 9 

+ ^ [lp' sin F{G - a + 2 - 2cosF) + ip sin F{S,[ - G') 

+<pF' (2 + 2 sin^ + (a - G - 2) cos F - 2(^1 + 6))] } • (42) 

Again, by variation the profile functions 6j 6 8'iid ip are obtained. The classical profile functions 
serve as source fields. 

The quantities and x are additionally functionals of the heavy meson fields. For the bound state 
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in the P-wave channel one obtains upon substitution of the ansatz H17(l 



drr'^ $'2 



2" 



:Rc 



$2 _^ 





2 1 








^ 2 



Ra 

r 



*5 



+ - 



-sini^ 

F' - 
H — 
r 

2V2cM 



Ra 



*2*o *o*i - e - 77^ 



gmv 
4V2c f 



^ (e - *2 - (1 - cosF) ^-0*2 

2C" C 
2tj'*o*i *i*2 + (G + 2) *2 



2Md 



1 



a 



[G (G + 2) *i - rG'^'2] $ - 



Rg 
2 



sinF 



^2 



[G (G + 2) + G' (2 + R^) $] *o 



(43) 



Here a prime indicates a derivative with respect to the radial coordinate r. Furthermore the abbrevi- 
ation Ra = cosi^— l + a{l + G — cosF) has again been used. The functional 1^ leads to the equations 
of motion for the profile functions <&, ^Poi ^ij and 5*2 of the fluctuating heavy meson fields. In these 
equations the classical fields generate the binding potential. The solution to these equations pro- 
vides the bound state energy, e and the bound state wavefunctions that are subsequently normalized 
according to eq. ifT^ll . 

Finally, we present the explicit expressions for the hyperfine splitting parameter for the bound 
state in the P-wave, c/. section 4. For convenience we employ additional abbreviations with regard 
to the light meson profiles defined in eqs JHJ and (|22|l 

Vi = cosF - a (Ci - 1 + cosF) , 
V2 = ■ 



The explicit expression for the P-wave hyperfine parameter, which enters the mass formula for the 



16 



even parity heavy baryon 131(1 . reads 
2 



X 



3 Jo 



dr r Px{r) 



2a 



e--cj) iV2-2Vi)-^i2 + R^)^ 



+ {2Vi + V2) [(e - ^-1 - % 



1/2*2 H </?*o 

r 



1 

~2 

H ip-^i % 



1 



2-^ *0 

r 



"2 -^2 i 

r r 



(2 + + AMd sini^<I>^'o 



d 
r 

Ay/2cM 



<! sinF 

r 



(2 + i?Q + Vi) 1'i*2 ¥'*o1'i 



-1/2*2 + 2a(^vI/Q\p2 



4V2c J" 



|(3Ci + Ca) *o*i + - (2 - 2^1 - 6) *o*2 + -^'*i*2 + ^V*2 j 
r r ) 



a 



.^^j + * + h/i - 



1^2 



c c 

(G + 2)vI.i_ 



(44) 



*5 



(2 + (2a +6-2) $*i + ^ [2^$' + (2 + Ro) - aGV*] *o 



(G + 2) 6$' + ( 1 + ) (ei + ^2) - c.^'v^ 



*5 



(45) 



Substituting the bound state profiles as weU as the soHton yields numerical results which are used to 
compute the heavy baryon spectrum according to eq. ((31|l . 
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